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ABS IKACT 


The optimum design problem of R.C. Caimney is 
formulated as a multistage optimization problem. The 
resulting non-linear programming problem is solved by 
Sequential Unconstrained Minimization Technique (SUlffll) 
using Quadratic Extended Interior Penalty Punction 
Method. Modified Uewton's method is used for unconstrained 
minimization. 

Extended interior penalty function technique for 
converting constrained optimization problem into sequential 
unconstrained optimization problem and Modified Uewton’s 
method for solving unconstrained minimization problem are 
discussed in detail. The analysis and design procedures 
for R.O, Chimney are briefly described as, by now, these 
are well documented. 

The efficiency of the methods used in the present 
work is established by comparing the results of and 
the C.P.U, time required to solve the same problem as 
given in Ref . (2, 11) where different techniques have 
been adopted to obtain the optimum solution. 



O^^TER 1 


nriRODUcrioH 


1 ,1 Introduction: 

Chimneys are relatively tall structures subjected 
to three t 3 rpes o±' stresses: (i) stresses due to self 
weight (ii) stresses due to wind or earthquake load and 
(iii) stresses due to temperatuce variation between the 
inside and outside of the chimney. Brick chimneys are 
suitable only for short heights, as they become bulky 
with increase in height and require heavy foundations. 

Also due to large temperature gradient, brick chimney 
frequently cracks, and becomes unstable. In contrast, 
concrete chimneys are lighter and stronger, and are less 
vulnerable to cracks due to temperature difference. If 
the temperature of the flue gases does not exceed 400°C, 
concrete chimneys can be used without any special fire 
brick lining, for higher temperatures, fire brick lining 
is provided with an air gap between the inner face of 
the chimney and the lining. 

Most reinforced concrete chimneys have a circular 
cross section for the reason that a circular section 
affords a minimum of friction which reduces the draft. 
Moreover, the circular cross section is easy to construct. 



2 


i'he materials used in the construction of a concrete 
chimnev are concrete, reinforcing ste-'-'l and refractory 
bricks. Any saving in the quanGity of materials used 
would result in the reduction of cost of chimney. An 
optimum design procedure for reinforced concrete circular 
chimneys is presented in this work with the goal to get 
a minimum cost design. 

In general the design of reinforced concrete 
chimney is a trial and error process. Normally the 
chimney is designed by assuraing mean diameter of the 
shell, thickness of the shell and percentage reinforccjment 
a various sections and then the stress analysivS is carried 
out to check if the stresses due to dead load, wind or 
earthquake load and temperature variation are within 
the permissible liraits. lenerallj^ such designs are 
conservative inyol’'.''ing repetitive computations and quite 
time consuming. I'herefore, it is natural to carry out 
computcsr aided design of such structures. As already 
mentioned, the ra?-'iin forces acting at a cross-section of 
chimney are dead load, wind or earthquake load and comperature 
effect. fhe first three t3rpcs of forces depend directly 
on the dimensions of the chimney above the section, therefore, 
a reduction in the size of the structure would in turn 
reduce the forces acting on it, Ihus the problem is most 
suitable to be formulated and solved as an optimum design 
problem. 



..he op timum design o l a chimney in the present 
work is cast Ci-s a mathematical prO'-gramming problora. Ihe 
objective is to minimize the material cost of the super- 
structure subject to the behaviour 'and side constraints. 
The possible design eniriables are the mean diameter of 
chimney shell, thickness of the shell a,nd reinforcement 
at various sections along the height of the chimney. 

As per code provisions six stress constraints are to 
be satisfied at each section, the number of design 
variables and the number of constraints increases directly 
with the number of sections to be checked and thus the 
size of the problem becomes large. 

Chapter 2 describes the .problem formulation. The 
optimization technique to solve the problem is described 
in Onapter 3. Results of a 80 m tall R.C. chimney are 
discussed in Chapter 4 and conclusions drav^n from the 

present work are also summarized therein. 

^ Rrevlous V^ork : 

The forces acting at a certain section of the 
chimney depend only on the geometry of the chimney above 
the Section. This motivated Kapoor and liariharan |2j to 
decompose the optimum design problem into a number of 
small sized optimization problems. In their formulation 
the chimney was divided into certain number of segments. 



starting from the top section. At each segment a nonlinear 
programming problem v\ras optimized via SUM! using interior 
penal fimcL'ion approach and Davidon 'fletcher Powell 
method for unconstrained minimization. 

Subsequently Bandyopadhyay |lll handled the 
optimum design problem of R.O. chimney via djmamic 
programming and showed sufficient saving by way of O.P.U. 
time, thus establishing the superiority of the d 3 mamic 
programming approach over SUMT for seeking the optimimi 
solution of this class -of problems, 

1 . 3 Present Work; 

I’he present work addresses itself 

(i) to solve the multistage optimization problem of 
R.G. chimney by converting the resulting constrained 
nonlinear programming problem into unconstrained 
minimization problem through extended interior 
penalty function approach and solve the resulting 
unconstrained problem by Modified Newton’s Method ; 

(ii) to compare the optimtim design so obtained with 
the optimum solution reported in Ref. \2, 11 j and • 

(ill) to compare the relative efficiency of the 

optimization techniques used herein with the other 
two techniques used in Ref. . |2,1l| for this 


class of structure 



D 


C'tie cost of concrete and steel are only considered 
in the objective function as in the earlier two works cited. 
Furthermore, the Indian Standards Specifications ll,3l 
have been follow^ed in formulating the constraints so as to 
compare the results of the present work with those reported 
earlier. 



CHAPTER 2 


PROBLEM PORMTJLATIOH 

2 . 1 Pormulation of the Problem: 

'?he optimum design of a chimne:/, in the present 
work, consists in minimizing the material cost of its 
supers bructure subject to the constraints imposed by the 
code provisions jl,3|. The resulting mathematical 
prograiiiming problem is to find the design vector X 
which minimizes the objective function 

C = f (X) 

subject to behaviour constraints 

g^ (X) > 0, D = 1 , m (2.1) 

and side constraints 

1 . ' u . . 

X 4 < X . < X . ,1=1, n 

The objective fimction is taken as the cost of 
concrete shell and vertical reinforcement used in the 
superstructure. The mean diameter, thickness of concrete 
shell and percentage of vertical reinforcement are taken 
as the three continuously varying design variables. In 
terms of the varis-bles , the cost of concrete shell can be 
expressed as 

H . 

C = / C n d(h) t (h) (1 + Sp(h)) dh (2.2) 

o . 



where d(h), t(h) and p(h) are the mean diameter, 
the shell thickness and percentage of reinforcemenc at a 
section of distance h from top of the c limney. K is the 
total height of the superstructure, C-^ is tne cost of 
unit rolume of concrete and S is the ratio of the .jost of 
unit volume of steel to that of concrete, 

Bj converting the integral in Eq. (2.2) to a 
summation form, we get 

ir+1 

0^ = 0^ Z m dj_ t^(l + S p^) h^ (2.3) 

i=1 

v'/here the index i denotes the ith section of the 
chimney which is divided along its height into W equal 
parts. It is assumed that the mean diameter and shell 
thickness of the chimney vary linearly between two 
consecuti/e sections. 

Sq, (2.3) can be rovvricten as 
E+1 

0^ = 2 hi , (2.4) 

i=1 

where Aj_ = tc dj^ t^(l + S p^) 

and hj_ is the disuance beuv/een the ith and (i-l)th 
sections (Eig. 1). Eor i=1 i.e. the top seccion, the 
mean diamete-r is fixed from nonstructural design 
considerations, fhe thickness of the shell, t^, and 
percentage of steel, p^i , at the top section have been 




FIG'1 THE CHIMNEY DESIGN PROBLEM DECOMPOSED 
TO A MULTISTAGE DESIGN PROBLEM 


set. to the minimr'm. The resiilting optimization problem 
-l.I “:-,o7e 3-1'^ design variables, viz,, and .Pj_, 

i=2, K+1. As the size of the problem tarns out 

to be large a kind of decomposition procedure was adopted 
in Ref. |2j and the same is being followed in the present 
v/ork. The dead load and wind load at any ith section 
depend evidently upon the dimensions of the chimney above 
that section. Hence it can be argued that if one proceeds 
to optimize the equivalent area of cross section, A^, 
sequentially, starting from the top, the I'esulting design, 
though not the theoretical optimum design of the entire 

structure, still approaches it from the upper side. 

H+1 N+1 

Min 0^ = Min (C^ 2 ^i^i) < ^ hj_(Min (2.5) 

i=1 i7=1 

Chis formulation consists of sequentially 
minimizing A^ = 2, H+1 , subject to the behaviour 

constraints on the chimney and side constraints on the 
design variables, for each section, the objective 
function is the equivalent area of cross section of 
concrete 

Af = •n; d^ t^ (1 + S p^) (2.6) 

A^ represents the equivalent area of concrete 
used in the computation of the cost of the superstructure. 

The problem has thus been decomposed to N problems 
each with three design variables dq_, t^ and Pj_» Tbe 



behaviour constraints are the bounds on stresses in each 
o-P the decomposed problem and the side constraints are 
the bounds on design variables as specified in Indian 
Standard Specifications |l|. 

2 • 2 Code Specifications; 

2.2.1 Loading conditions : 

The Indian Standard Specification for the design 
of R.O. Tnimney jl| requires that the chimney section at 
all heights should be safe against the following load 
conditions. 

(i) Dead load + v;ind load 

(ii) Dead load +• earthquake force 

(iii) Dead load + temperature effect 

(iv) Dead load + v/ind load + temperature effect 

(v) Dead load + earthquake force + temperature effect. 

2.2.2 Permissible stresses; 

The permissible stresses in concrete and steel 
for each of the above loading cases are specified in the 
code 1 1 I and are given in Table 1 for the sake of 
completeness of presentation. 



.']A.BLE 1 


PERMISSIBLE STRESSES 


Load condition Permissible stresses 

(kg/cm^) 


Concrete Steel 


(i) 

Bead load +wind load 

0.38 

^c 

0.57 

^sy 

(ii) 

Dead load + earthquake 
force 

o 

• 

o 

^c 

0.60 

fsy 

(iii) 

Dead load + temperature 
effect 

0.35 

fc 

0.55 

fgy 

(iv) 

Dea-d load + wind load 
+ temperature effect 

0.50 

fc 

0.65 

^sy 

(v) 

Dead load + earthquake 
force + temperature 
effect 

0.30 

^c 

0.65 

fsy 


where is the 28- days cube strength of concrete in kg/ 
fgy. is the yield stress of steel in kg/cm^. 



2.2.3 I/Iinimum requirements; 

(i) wiinimum thickness of vhe shell (tj_) = 15 cm 
' for -diameter (d^^) _< 6 m 

di -6 

t > >15 + — r-r cm if d,- > 6 m 

^ ~ t,2 ^ 

(ii) Percentage of reinforcement 

0.003 < Pj_ < 0 .05 

(iii) Minimum clear cover = 5 cm . 

2.2.4 Stress calculations : 

At any height of the chimney the dead load, wind 
and earthquake load can be calculated once the geometry 
of the ciiimney above the section is known and the stress 
in each case can be computed as discussed in the succeeding 
paragraphs. Among the wind load and earthquake load, only 
the more severe of the two a cross-section* is considered 
for ;hG purposes of stress analysis. 


* It was observed by Kapoor and Hariharan j2|, that the 
wind loads are in general predominant in the lower 
half of the chimney while- the earthquake loads dominate 
in the upper half. Ihe nature of wind and earthquake 
moments along the height of the chimney considered in 
the present' work is sho^ivn in Pig, (2) which is taken 
from Ref. (2) , 





(i) Stresses due to dead load only 


the maximu’Ti vertical stress in concrete due to 
dead load only is given by. 

(a) For annular section with no opening 

W 


f ’ 

^ cvd 


2ixrt 


(2.7) 


(b) For annular section with one opening (Fig. 3) 

W 


■^cvd 


( 2 , 8 ) 

2(7i-P )tr 

Ihe maximm vertical tensile stress in steel is 


given by, 


■svd 


m f 


cvd 


(2.9) 


( i i ) Vertical stresses due to temperature only 


Ihe difference in temperature betv/een the two 
faces of chimney shell is gi'ven by 


t D- 


T 


■ n 


bi 


X 


Cp Dp 


1 


+ — 


^b ^bi 


Di 


'bi 


^%i ^ D]3j[_ 


! r^k-1 r^^ Ci^ D^ ^3 Dg ^q-^c ^2-^co 


G koD. 


"q b b 

( 2 . 10 ) 

where f is the maximum temperature of gas inside chimney 
Tq is the minimum temperature of outside air • 
surrounding the chimney. 

^q» ^c» *^b’ » ^2 ^s 2i.re coefficients depending 



on thermal properties o-^ concrete, lining and insulating 
medium (air). The geometric parameters appearing in 
Bq. (2.10) are shown in Big. (3). 

The vertical stress in steel, fgijiy » concrete, 

foTY’ temperat;n?e gradient through the shell 

(Big, 5), are given by 




= a (3 - K ) Eg 

(2.11) 


fQIJY 

= a K Eg 

(2.12) 

and 

K 

= pm + Y'pm(pm+2^) 

(2.13) 

where , 





a is the coefficient of thermal expansion of concrete 
and steel per degree centigrade. 

3 is the ratio of distance between inner surface of 
chimney shell an^l vertical reinforcement to total 
shell thickness t. 

Eg is modulus of elasticity of steel, and 
Eq is modulus of elasticity of concrete 
m is the modular ratio 

p is the ratio of total area of vertical reinforcement 
to total area of concrete chimney shell at the 
section under consideration, 

(iii) Moment due to earthquakes 

fhe procedure for calculating moments due to earthquake 
forces for stacklike structures is laid down in section 5, 3 
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of earthquake code |3| , 


The approximate fundamental period of free vibration, 
T in sec. is calculated from 


W H 

T = Co, Y (2.14) 

^ E A-g 


where , 


C,j = A coefficient based on slenderness ratio of 
chimney given in Ref. jjj. 

W = Total weight of the chimney including weight of 
lining above the base. 

H = Height of chimney above the base 

E = Modulus of elasticity of material of the 
structural shell 

= Area of cross-section of structural shell of 
the base 

g = Acceleration due to gravity. 


Using period P, the horizontal seismic coefficient 
is obtained from 


a 


h 


where, I 

P 


-"0 


Sa 




(2.15) 


Importance factor for the structirre • : 

A coefficient depending on soil foundation system; 


Basic horizontal seismic coefficient given 


in Ref. \ 3 \ | 

Average acceleration coefficient from average 
acceleration spectra given in code and reproduced 
in Eig. (5) • 




-Phen the design bending moment iV^ at a distance x 


from cop is calculated by the followin'? formula 




0.6(x/H)^/^+0 



(2.16) 


where 

h is the centre of gravity of structure above the base, 

( i V ) Stresses due to lateral load and dead load : 

-'he whole section of shell is under compression, if 
(a) . for anniiLar sections 


if e/r < 1/2 (2.1?) 


The maximum vertical compressive stress in concrete 
shell is given by 


f 


cw 




27irt 




(2.18) 


where, 

e is the eccentricity 

r is the mean radius of the shell at the section under 
consideration 

t is the thickness of, the shell at the section under 


consideration 

W is the weight of the chimney above the section, 
(b) For annular section with one opening (Fig. 4-) 


e 


1 


< 


2 ) 


(tc-P )^-sin^^ 
(7C-p)co3 P+sinp 


\ 

“ 3 sin§ I (2,19) 



iL\J 


I'lien in sucLi case, the maximum vertical compressive 


stress is concrete is given by, 




■ cw 


2 (tc-P )r t 


1 + 


r e sinp 7 f ■> 

^ j )cos|3+sinP I 

(ic-P) - 0.5 sin2p - 2sin^/(K*.-P ) 

( 2 . 20 ) 


If e/r is greater tnan the corresponding right 
hand side of expressions (2,17) or (2.19), tbe angle a (Fig. 4 ) 
defining the position of the neu'cral axis of the R.C. shell 
section is calculated by numerically solving the transcendental 
equation given below 
e k 


r 



( 2 . 21 ) 


where , 

A = (l-p)(a-sina cosa)-(l-p+mp)(P+sinp cosp 

-2cos^sinP) + mpTC 

B = (l-p) (sina-acosa)-(l_p+mp) (sinl3-Pcosa)-mpTECOsa 


m is the modular ratio 

l3 is angle in radians of semi-flue opening 
p is ratio of total area of vertical reinforcement 
to total area of concrete chimney shell at the 
section under consideration, 

Regula; falsi method has been used to find the value 
Once a is Imown, the compfessive stress inconcrete 


of a . 



at raid fibre 


•f ’ 

- cw 


W 


of the section (Pig. 4), is calculated as, 

(cosp-cosa) 

(1-P)(s ina-acosa) -(1 -p+mp) (sinjS-pcosa) -mpitcosa 


( 2 . 22 ) 


The maximum compressive stress in concrete, fQ^f becomes, 


f 


cw 


f ' 

- cw 


1 + 


t 


2r cos|3 (cosP-cosa) 


(2.25) 


L’he maximum tensile stress in steel, fg^y> is given by. 


■f = rnf* ^ 


1 + cos a 


cosp ~ cos a 


(2.24) 


(v) Stress due to dead load and temperature 

In the case of reinforced concrete, the stresses 
cannot generally be added arithmetically but must be 
combined according to the elastic properties of the materials 
and in the proportion of reinforcement to concrete. 

The maximum vertical compressive stress in concrete 
due to combined effect of vertical load plus temperature 
effect is given by, 

^cvdt “ ^CTV 


K 


if %t < 1 


(2.25) 



0 r 


f 


f, 


;Vdt 




OTY 

K 


2nip3 + 1 


2 (1+iap) 


i±' K^i- > 1 


(2.26) 


where, | 

= -mp + Ymp(rap+2Z)+2K(l+mp) 


^cvd 

-CTV 


is the ratio of distance between inner surface of the 
chimney shell and the neutral surface resulting from , 
combined dead load and temperatiire , to the total shell 
thickness, .t. 

The maximum vertical tensile stress in steel due to 
combined effect of dead load and temperature is given by, 


f 


^svdt 


STY 


(Z-K) 


Z+mp - (mp (mp +2Z ) -2mp ( 2 -K ) 


fsvd 


fgrjiY 


) 


1/2 


(2.27) 


( V i ) Stresses due to dead load, lateral load and temperature : 

L?he stresses in the concrete and reinforcement can 
be computed for each case, i.e. on windward side and leeward 
side of the cnimney. It is clear that the desi,gn will be 
governed by the worst cases, namely (l) compression in 
concrete on the leeward side of the shell and (2) tension 
in the reinforcement on the windvjard side of the shell. 

The maximtim compressive stress in concrete due to 
combined effect of lateral load, dead load and temperature 
figs. (6 and 7 ) is > gi'ven by 






Stress 'iue to 
D L and .ateifd 
load 


Stress due to 
temperature 


X ■. 

^sv-comt» 
T m 



lateral load 


DN THE LEEWARO 




! Stress due to ,, ^ ‘ 

^ weight and ' i I 

laterc'i load I 
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Stress due to lateral load, 
weight and temper 0 tufe 

FI0.7 STRESS DISTRiBUTlON ON THE lEEWARO 
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iE 

■U 


m 



o 
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St f ess due to 
lateral load 
and weight 


Stress due to 
temperature 



Stress due to 



load, weight and 


temperature 


FlG-8 STRESS DISTRIBUTION ON THE V/INOWARO 
SIDE OF THE CHIMNEY 



f 


cw-coinb 


or 


4 ^ 

cw-comb 


■CJY -^oomb 
K 


^comb — ”* 


(2.28) 


1.Q2Y 

^ T 

OW ^ |r 


2mp2 + 1 

2 (1 + mp) 


^‘^comb — ^ 


(2.29) 


'^'/here is the ratio of distance between inner 

Surface of chimney shell and neutral sur-^'ace resulting from 
combined lateral load, dead load and temperature effect to 
shell thickness, and is calculated as 


^C'omb ~ "P™- I (pm+22) + 2K (1+pm) 


f ’ 

^ cw 


^CTV 


1/2 

I 

(2.30) 


fhe maximum tensile stress in reinforcing steel, 


^sw-comb ’ windward side of the chimney, (Fig. 8), is 

given by 


-P 

^ s w-comb 


7z-K) 


( 3 +pm) 


(pm(pm+2Z) 


^sw 

-2pm(Z-K) — — 
I STY 

(2.31) 


1/2”^ 

) 


2,3 Compatibility Requirement; 

For the purposes of a compatible and practical 
structure, the mean diameter plus the thickness of the 
shell at ith section should be greater than or equal to 
the mean diameter and the thickness of the shell at (i-l)th 


section i.e., 



2 . ^ Optiim.m Design Problem: 


-lie minimum cost design of R.O. chimney has thus 
been formulated as a multistage optimum design problem. 

It consists of finding the solution of IT optimum design 
problems each consisting of 3 design variables, 6 behaviour 
constraints and 4- side constraints as follows: 

minimize Aj_ ( X ) = -n: dj_ t^ (1 + S p^) 

i = 2, IT + 1 

where X = (d, t, p) 

Subject to behaviour constraints ^ given by the 
limitations on stress defined in sections, 2. 2. 2 and. 2,.2,4 
and side constraints defined in sections 2 . 2.3 snd 2 , 3 , 

2,5 ITormalization of Constraints ; 

All the constraint equations are normalized such 
that they are non-dimensional and lie in ilie range of 0 to 1 , 
Typical examples of normalizing a behaviour and side 
constraints are shown below. 

A behaviour constraint is typically represented as 
S < ^max 

Rewriting the above inequality 


^max “ ^ ® 




the side constraints also. In the present formulation 
all the constraints (g-;S) are expressed in greater than or 

cl 

equal to zero form. 
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5 • Introdac slon: 

In the interior penalty function formulation for 
solving inequality constrained minimization problems, the 
constrained minimization problem is transformed into a 
sequence of unconstrained minimization problems. It is 
reliable and it facilitates the use of efficient unconstrained 
minimization algorithms to solve the inequality constrained 
problem. Moreover, in structural design optimisation the 
method exhibits a desirable property, from a practical 
viewpoint, of generating a sequence of feasible designs 
that bend to funnel down to the middle of the feasible 
region. As a result of this property each design generated 
defines a structure for which the excess weight, with respect 
to the minimum, is distributed so that none of the constraints 
become critical prior to converging to an optimima design, 

■®e interior penalty function is only defined 
within the feasible domain. It is known that if the starting 
point is feasible, then the method should generate only 
feasible points in the sequence of desi,gns converging toward 
the optimum. However, most of the powerful algorithms lor 



unconstrained minimization involve a sequence of one- 
dimens lonal minimizations and in practice they are carried 
out oy the numerical methods which may well generate trial 
desi m ]points outside of the feasible region. Extending 
the penalty, function formulation into the infeasible 
design space provides an efficient mechanism to recover 
from violations brought about by analysis and also permit 
the use of an infeasible initial design as a starting 

# 

point. Numerous recovery schemes that can bring a design 
from infeasible design space back to feasible design space 
are available. But it is believed, Ref. 114| , that these 
recovery schemes are not as efficient as the extended 
penalty function approach. 

5 2 Motivation for Extended Interior Penalty Function Method: 

Oaroll proposed the interior penalty -function 
method go convert the constrained non-linear optimization 
problem into unconstrained problem. In general the 
constrained optimization problem may be posed as 

Minimize 1? ( X ) 

Subject CO constraints g^ ( X ) ^ 0 > ' j • • •^con * 

In the interior penalty function formulation a 
new function, penalty term is chosen such that , its value 

P(1 * r) i® toy augment ing a panalty t®ie» to 


til® otojooti’wi fnaotioa# 



will be small at points aijvay fpora th.e constraint; boundaries 
and will tend to infinity as ‘che constraint: boundaries are 
approached. Phe penalty function P ( X , r) is given 
by daroll 1_1 9 H 



■^'^con 

? ( X ) + r 1 ( X ) (3.2) 


.3(1) = 


1 

gj ( X ) 


for g^ ( X ) > 0 


Phe term r ‘jj(x) represents the penalty associated 
with the j zh constrain, and is an interior penalty function 
in the sense that it is defined only if X is inside the 

feasible design domain. With X denoting the point in the 
design space, P (X , r) attains its minimum value for a 
given value of r. It is shown (Ref. 8) that as r goes to 
zero 

min. P (X,r) - P (T* ) 

3*ncl X -*■ X • 


■Oherofore, the penalty function is sequentially 
minimized for a decreasing sequence of r until r approaches 
zero, 'Phis is called sequential unconstrained minimization 
technique (SUIT) and uses any of the known unconstrained 


minimizers . 



:?he variation of the ordinary interior penalty 
function along a direction S in the design variable space is 
shovm in tig. 9, in which a represents the scalar variable 
that denotes the distance along S . It is assumed that the 
minimum point M in the direction S ,• is to be found by a 
numerical method, say, cubic interpolation. It will 
therefore, be necessary to locate two trial designs one on 
either side of uhe minimum M', Steps are to be taken in 
the direction of increasing a to find a point on the portion 
Mk of the curve. Ihis search becomes an increasingly 
difficult problem because the interval M'N' tends to become 
very narrow as the SUB/ff algorithm advances toward the 
optimum design, fherefore, it is likely that the search for 
a values in the interval 'hV N* will produce trial design 
points that are outside of the feasible region. Iechniq.ues 
for stepping back into the region M' N’ usually require a 
large number of function evaluation. Every function evaluation 
entails structural analysis and thus the procedure of stepping 
back introduces a heavy computational burden. The troubles 
do not end once a point in the M'h' region is found. This 
is because the function value and the first derivative of 
P(X,r) with respect to a will frequently be very large at 
the located point in the interval M'!’.. These high values 
are not well suited to construct a good cubic approximation 
of P (X,r). Therefore, the minimum of the cubic approximation 
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PIG 9 ONE OIMENSIONAL MINIMIZATION OF THE 
INTERIOR PENALTY FUNCTION 



does not yield a good estimate of the a value corresponding 
CO the actual solucion of bne one-dimensional minimization 
problem. 

fhese important drawbacks motivated Oasis and 
Schmit 1 Ref. 13 ] , the idea of devising an extended 
penalty function which is defined in the infeasible region, 
while still retaining the advantages of an interior penalty 
funccion formulation. The extended penalty function 
formulation accomodates the infeasible design points, 
avoiding bhe costly search for a feasible design in the 
narrow region M' R'. 


3 , 3 Linear Extended Interior Penalty function; 

Kavlie and i4)e j__ 20 ^ proposed a formula 

for the extended penalty function. The transition point 
in the interior penalty function satisfies continuity of 
both function value and its first derivative. They define, 
the individual constraint contribution to the penalty terra 
of the compound interior-extended penalty function as 
follows 



(3.3) 



where is a transicion point between the two 
con>s "’raiiic; functions and varies as 

So = Cr (3.4) 

where 0 is a constant and r is the penalty parameter. 

■Because the G-^ of Eq . (3.5) have discontinuous 
second derivatives at the transition point, the linear 
extension of :;he penalty function is not suitable for a 
second order (e.g. hewton’s method) optimization algorithms. 

fhe formula proposed by Kavlie and Moe (Eq, (5.5)) 
defines the extended penalty part of each constraint gj. 
i'he ke:/ parameter gg corresponds to the values of the g!® 
at which the transition between the regular inferior penalty 
function definition and the extended penalty function occurs, 
From the one-dimensional view of the interior penalty function 
shown in Pig, (9), it can be seen that definition of a 
useful extended penalty fiinction requires the transition 
point to lie on the Mk portion of the curve, the reason for 
this requirement is that the extended interior penalty 
function would not include the minimum at point M and indeed 
its minimum (Mg) might even be outside of the feasible 
region, as shown in Eig. (10). Ifow, recalling that the 
introduction of the extended penal cy function was motivated 
to locate the minimum along linear direction, it is obvious 



that extended penalty function should be defined so as 
to include the minimum of bri'^inai penalty function, A 
more rational approach to find is described in (Ref.13)» 
Hov'/ever, since it is not pertinent to the present work, 
it is not being discussed here* 


5 • 4 Quadratic Extende d Interior R-enalty -Ttinction: 


The disadvantage of the discontinuous second 
derivatives ■■at the transition point in the linear extended 
penalty function can be overcome by introducing a quadratic 
extended interior penalty function (Ref. I 4 ). It is 
continuous and lias continuous first and second derivatives. 
'The definition of in Eq. (3.2) for a quadratic 

extended penalcy function is given by, . ■ 


1 


if gj ) >: go 



£ 2 ® 

go 


2 . 

1 ) 


So 


- 1 ) +1 


(3.5) 


fhe transition point for the above quadratic 
extension to ’the penalty function is assumed (Ref, ) as 

g = G r^ ; 1/3 :< P < V2 (3.6) 


where C is e* consta-iit 
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As r goes to zero, the above equation has co satisfy 
the conditions: 

(i) r dj(x) -* Q for any >0 » 

which represents a vanishing contribucion in the 
feasible domain and 

(ii) r G-^(X) -> CO for any (X) < 0 , 

V'f'iich represents an increasing penalty for violated 
constraints. 

the lower bound on the exponent in Eq, 3.6 guarantees 
that i:he penalty for violated constraints increases as r 
goes to zero while the upperbound on the exponent is needed 
to keep the minimum point in the quadratic range of the 
penalty function, the quadratic and linear extensions to 
the penalty function are shown in Eig, 11 along with a 
conventional interior penalty function for comparison. 
representation of the quadratic extended penalty fuaction 
for different values of r is shown in fig. 12. 

3 . 5 Avoiding Negative Design Variables: 

It is to re-emphasize that the extended penalty 
function resolves the difficulties arising from the fact 
that the. ordinary interior penalty function is defined only 
in the feasible region. Even though extended interior 
penal -cy ^unction is defined in the feasible region, it is 



Constraint function 
FiG^II PENALTY FUNCTIONS 


^ >r2>*’3 



W' I 

I I 

I I 


I I ' 

» I ' 




valid only over that portion of the design space v/here 
tae :.;onB train :s are meaningful. For example, g!®loose 

t] 

their pir/aical significance Vilien one or more of the design 
variaoles calces on a negative value. Therefore if a search 
for tiio one -dimensional minimum leads i:o points wirh 
negative design variables then analysis should not be 
por.iorraed there. .Che search instead be returned to the 
porcion of the design space where all the variables are 
posicive. The method in jRef. 13 j , resolves this problem 
in one step, eliminating the uime consuming task of making 
several trials. :f'or completeness of presentation it is 
described h-erein. 

Fig. 13 shows a one dimensional view of extended 
interior penalty function. Che curve represents a plot 
of P for designs located along a straight line passing 
tnrough the point and having direction S . For illustrative 
purposes the design space has been assuriied two-dimensional, 
as s.'io’m in Fig. I4. Chere, the location of straight line 
in the direction S can be seen. Several points have been, 
marked on this line in .Figures I3 and H, their meaning is 
described in the next paragraph, 

Ji-ssume that starting from a point A in one dimensional 
minimisation in the direction S , determined by a certain 
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is to be performed. -Phe starting point A is 
kno .11 and ii is assumed to be the minimum point along 
the previous direction, Ihe scalar variable a locates design 
poinis a, long the direction B. Phe scalar 'a' has its own 
origin au point ’A' and ir is taken positive in the direction 
of decreasing P at point 'A 'along S . Phe minimum point 
sought along 3 is designated as point M, At point 

B along the direction S the constraint gj^ = 6 and therefore 
it defines the va.lue of a corresponding to the transition 
point of extended interior penalty function. At point 'C* 
the constraint gj_ = 0, therefore, this point is on the 
boundary separating the feasible region from the infeasible 
region, at point 'E-j ' the lower bound side constraint 
corresponding to the design variable .x-i is equal to e and 
B 2 has similar meaning with respect to X 2 . At point E 
the side constraint related to equals to zero. Point 'G' 
where x-^ = 0, is the first point along the line for which 
I' is no.; defined, i^t points X and Y the function P is 
undefined. Point X is characterized by the fact that 
x^ < 0, X 2 > 0 while Y is located in the region where x^ 
and X 2 are both negative. 

Starting from A, the search seeks a point on the 
other side of M such that an interpolation function can 
be constructed and an approximate’ minim’um may be determined 
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from io. If step si'^e is such thac either X or Y is 
Picheci ur , the analysis used to test whether or not the 
poin: is usetul '’or the interpolation cannot be performed 
and another point has to be tried. Jor the present purpose, 
the poinos 1, ? , .1-^, 0 and B, located on opposite side of 
ivi with respect to A are prospective candidates. Of these, 
point 1 can be excluded at the outset since ic corresponds 
to x-| = 0 and deletion of eli.mencs is not considered by 
current analysis. ..’'urt heritor e, points B and 0 are difficult 
to locate because is usually a complicated implicit 
function. iluorefore, only points E-^ and E remain viable 


candidates. 

It is straight forward task to locate either point 
E-] or E because they represent the intersection of a 
straight line in direction S with hyperplane. However, B-j , 
has the advantage of being closen to M and therefore it is 
a better point at which to generate data for determining . 
an approximate minimum. -ho detailed procedure for locating 
B-j is now described, a lower bound side constraint function, 
related to a design variable Xy is witren as 





(3.7a) 


where 


represents the lower bound allowable value 


of the variable Xj. 
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iit point 51^, = e, so the 3 th co-ordinace 0 ! die 

point 3j can be evaluated from Equation ( 5 . 7 a). 

^^2 ^alj^^+e) (3.7b) 

low assume that starting from a feasible point 
a step size ay is applied in the direction of S to get 
a point Y having some negative components, let f denote 

the set of all subscripts corresponding to negative design 
variable components at point Y, and let 'm ' be number of them. 
Eirst, point Y is determined by 

x(Y) ^ y(A) + cc g (3.7c) 

laving identified all the design variables that 
are negative at Y, there will be 'm' points (j 6 3 ) 

to choose from, yet only one of them will satisfy the 

(E .) 

condition of belonging to the region where all the Xj_ 0 ,s 
are greater than zero. From Fig. ( 13 ) it is easy to see 
that the correct is the one nearest to point A. As 
this distance is measured from A by the value a^, the 
problem reduces to selecting 

“BK “gs “Ej e.7a) 

•)(* 

Af ter all (3S 3 ) have been computed. 
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.’o obcain the s consider the relation. 


f Ciil -t ) 

d 


= + a,, S 


(3.8a) 


rio'vvo /er, the ,j ch component of is known, 

a^ ^.CiOv/ii oy -iij nation 0. 7 b) for tnis component equation (3 -8a) 
imnlies 


= .f) sj 

and equasions (3.7b) and (3.8b). yield 


(3.8b) 


^31 


^alq (1+e) - 


(A) 


(3.8c) 


■it :;er f ini in 3 a;^ from the condition e.xpressed in equation 
( 3.Td) , the sou'dit a-^ter point is determined by 


X 


Y (■‘^ 4. ^ ^ 

A + S 


( 3 . 8 d) 

, ( 3 . 7 d ) and ( 3 .Sd ) show that it is 


Equations (3.8c) 
simple to amoid design where the analysis variables are 


not defined, 

5 ^ Initial 7 alue of r: 

.'ho initial value of the penalty parameter, r^, 

has to be carefully chosen for better convergence. It is 
customary to choose r^ such tiat the objective function 
and the pnmalty term are of the same order or magnitude 
at the starting point. Thus, 

I ^ 

^0 = I 2 i/g^ 
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5 . 7 ■ Gradie nt val uat ion.: 


jralienij '.or a quadratic e abended penalty function 

j. f:3 c c:i j r ■ , : 'i I: f,: a a, s 


51' 

a:^ 

r: 

neon g„ ^ 

ex. 

3x^ 

J- ^ 

3=1 ax^ 

!b = 

1 

agn(x) 

. Y w 


S^(x) 

ax. 


(3.10) 


> g, 


, 9gi(X) 

3 “ 5go) X — — 


gi(^0 < So 

(3.11 ) 


? X,- 


In the present work, partial deriyati^res , 
j are oaljulated by central difference mechod. 


3 . 8 .-iessjan Matrix Development': 

lessian matrix (H) is the matrix o^ second derivatives 


of penalty funct-ion, P, and it is given by 


p 

3 P 


( 3 . 12 ) 


a-nd 


? 

3 P 




ax^ axk 


neon 6 Ol 
+ r 2 


3X-:SX 


(5.13) 





‘'rom^ 3 qn. (3.5) , for tke quadratic extended penalty 


line : ) 1 , 


(j 


ox J 3X1^ 


a 


is given by 


cx c::j^ 


3 gl dg^ 

for 


a^gi 


ax^ax^ 


(3.14) 


s 


-3 ^ 

o I ax. axj, 


+ Sn( 


2Si 


- 3 ) 


a^Si 


for g^ < g^ (3.15) 

Because of the factors gT^and g“^in Sqs. (3.I4, 3.15) , 
the main contribution to the second derivatives of tne 
penalc'.'- I'unction is from the constraints which are nearly 
critical. For these constraints the following approximations 
for the second derivatives of the quadratic extended penalty 


function may be used, 


2 

a ^ n 

0 tTj 


dx_j 3x^ 


3 iSt 8 g, 

2 g. 

''='0 


= 2g; 


-3 9 St 9 gi 


axi, 


for gj. > gQ 

(3.17) 

for gj; < gQ 


It is noteworthy that this approximation requires only the 
first derivatives of the constraints. Therefore, the 
computational effort needed to generate the Hessian Matrix shall 



me tho d 


be much less. 

3.9. .bod i - led Jev/ton 's Me :hnb : 

3 • 9 . "1 Introduction: 

lie in .-one trained minimization od the penalty 
: unction, j:'(.(,r), is performed b^/ modi;''ied Hei-vton's 
in ulie present ’v<^ork. -Dhe disadvantage of the q_uasi— ITewton 
algoritcjns is that the number of function evaluations req_uired 
for the optimization procedure (which is measure of compu- 
tational ef ficiency of the algorithm) is a linear function 
O':’ the number of design variables. These algorithms, 
therefore, are not suitable for use waen a large number 
of desig-i variables are involved. It v\ras reported in 
Ref. 1 18} , that generalized Newton's method applied with 
an interior penalty function formulation can be used to 
overcome this disadvantage because it is possible to obtain 
a simple approximation to the second derivatives of the 
penalty "unction necessary for generalized Newtons method. 

This imitiiod r>,iquires a smaller number of function 
evaluations and is independent of the number of design 
variables. It is, therefore, useful for solving optimum 
structural design problems having a large number of design 
variables. However, in the present work, Newton's method 
with modification is used for unconstrained minimization 



with an intention ‘;o develop a general purpose computer 
P":’;) gramme .;o aanile large probleij^ as is common in structural 
analysis even thougi, due to decomposition principles 
applied, the problem in the present work does not have 
large number of design variables. 

There are two tneoretical arguments against the 
generalised liewton's method: 

(i) If riessian matrix (H) is not a positive definite 
matrix, a move in the direction given by, -{jfj 
mav result in an increase rather than a decrease 

in P (X,r), yielding, aj_=0, and terminate the process 
at the starting point Xj_. 

(ii) 'die Hessian matrix may not have an in verse, even if 
P(l,r) is convex. 

The modified second-order method used in the present 
work takes into account these limitations of che generalised 
Newton method. The direction vector S is generated according > 
to the algorithm given below. The basis for the algorithm 
is given in P.ef. jsj reproduced herein for the purposes 
of completeness of presentation. 

Che direction vector S is generated according to 
two rules. In both cases = ^i + “i ^ » where a^, is 

chosen co be the smallest value of a > 0 for which 



^ i‘- a, local minimum of ? f^r). 

■’he naes are as follows: 

'■" negative eigenvalue, let S be a vector where 

S" [h] S < 0 and Sv P < 0 (3.18) 

(ii) If H has all eigenvalues greater t nan or equal to zero 
caoose S such that 

either [h] S = 0 , S^yrP <0 (3.19) 

or [h] S = _yP (3.20) 

in which case, ifvT'P ^ 0, it can be proved that 

S^Q-I] S = -s Vp > 0 . Both Eqs. (3.19) and Eq. (3.20) 
can not hold. Let S satisfy (3.19) or (5.20). 

The onl7/- case in which Rules 1 and 2 fail to generate 
a non zero direction S is when [h] is a positive semidefinite 
mat3:*ix and "yP =0; that is, a point that satisfies the 
first-and second-order necessary conditions that be a 
local constrained minimum of P (x,r). 

:]he rationale 'for Rule 1 is that if 'the second 

\ 

derivative matrix has a negative eigenvalue there are certain 
directions along v/hich the function P(X, r) decreases and 
along which the rate of decrease also decreases. That is, 
for a vector S satisfying Eq. (3.18), 


and 


d P /da = ■S'vP <0 

d" I- (Xj_^^,r) / dx^ = 3^ [h] S < 0 (3.21) 

at a = 0 

''.'iius, unless ? varies along tne rav S emanating from 
so that e/encually d P (Xj^^^,r)/ da >, 0, the function 
value will tend to - This rule provides a guide for 
locating a region where Qf] is positive semidefinite, Por 
the minimization of a nonconvex function a few iterations 
along* vectors satisfying Eq. (3.18) can result in significant 
progress in minimization algorithm. 


The same remarks apply to (3.19) of Rule 2 , 


Mat ^ ^ 


0 holds for a indefintely large, 


P(Xi^ 1 ,r) will tend to - finally, if [hJ is positive 
de'i'ini'Ce , Eq. (3.20) yields the Eewton vector, when 
Eq. (3.20) is used, in general, 

S ■ = - 1h]"Vp ' (3.22) 


It is to be noted that any S' satisfying Eq. (3.21) 
has the desirable property that if P (X,r) is a positive 
semidefinite quadratic form having a finite unconstrained 
minimum. 


Rules 1 and 2, although quite useful computationally, 
are not guaranteed to generate a sequence of points having 


t 
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1 - 1 -nit poin-;s satisfying the first- and second order 
necofaSiary condi uions I'cr s- niniinuiii, A. cls-ss of direction 
vectors for v/dicii any existing linit point has this property 
is esuablisned in Ref. | 21 [ one particular e.xample of 
this type of direction vector is 

S = - A + 6^ e^ (3.23) 

where A is a positive definite matrix, e^ is an eigen 
vector of [^IlJ with minimum eigenvalue, and 6^ is a scalar, 

whore ^ 

I -1 if [h] has a negative eigenvalue and 

„i j P'-"' e^ > 0 

° 0 if [h] is positive semi-definite 

1 if [H] has a negative eigenvalue and 
V— "n - 

P e^ < 0 

-from the above observations, the suggesced algorithm 
is given as follows. 

5*9.2 Search Direction for Modified Newton’s Method: 

i'he searc.i direction S for one -dimensional minimization 
after developing Hessian matrix Qd is given by the following 
algorithm. 

(i) Reduce the Hessian matrix Qd to the- form 

H = L D L" ' (3.24) 

where 1 is a nonsingular lower triangular matrix with units on 
the main diagonal and D is a diagonal matrix. 





(iiy If D hias all posi-<:ive diagonal elements, solve for 

S = - ? (3.25) 

(iii) If D has some diagonal elements that are negative, 
solve L"" T = A . where A is a column vector with 
jth component = 0 if the gth diagonal element of 

D > 0 and with jth component = 1 if the jth diagonal 
element of D < 0 , Let S = T if _< 0, and 

S = - T otherwise. 

(iv) If 13 has all nonnegative diagonal elements, and- at 

« 

least one is zero, the vector S is generated according 
to S = - [h] ^ S/P . 

3 . 9.3 Algorithm for Step Leng:th for BJodified Mewton's Method : 

Quadratic interpolation method has been used to find 
the minimum of P (X, r) in a direction S . Ihis one- 
dimensional minimization method finds the minimizing step 
length a* ' in two stages. In the first stage, the S- vectpr 
is normalized so that a step length of a= 1 is acceptable. 

In the second stage, the function P (X,r) is approximated by 
a q_uadratic function h(a) and the minimum a* , of h(a) is 
found. If a* is not sufficiently closer to the minimum a*, 

the third stage is used. In this stage, a new quadratic 

2 

function h' (a) = a' + b 'a + c 'a is used to approximate 



fhis procedure 


i-'(\',r), and a new value of a* is found. 
i3 continued until a a* , which is sufficiently close 
to a* is found. fhe algorithm of this method is well- 
documented in text books and is discussed in detail in 



CHAPTER 4 


RESUT/T3 AND DISCUSSIONS 


4 . 1 Introduction ; 


An illustrative example using the techniques of 
optiirtum design described in previous, chapters is discussed 
herein. The problem, that was solved in Ref, [2 | , and 
later in Ref, j11 j has been taken up here too 
for the sake of comparison. The results are presented and 
discussed and conclusions are drawn. 

All computational work has been carried out in 
Tii-O 10S)0 system at I IT, Kanpur. 

4 . 2 The ^example Problem; 

The design of a 80 m high ventilated R.C. Chimney 
i3 cousidered here. The relevant data are given below: 

(1 ) Geometry: 

Height of chimney = 80 m 

Type - lined and ventilated 

Ventilated air space =10 cm 

Thickness of fire bx'ick lining = 11.5 cm 

The lining is supported at corbels at 8m 
intervals 

Hean diameter of shell at top = 4.85 m 
Shell thickness at ' top = 15 cm ■ 



Minimum vertical reinforcement percentage throught 
the chimney shell =0.3 

Flue opening of ■ 3m x 3m at 8 m from base. 
i'emTioratiire s 

Temp:: rature of flue gas = 160°C 
Mininram temperature of atmosphere = 4°C 


1 

= 0 

.5 ■ 


0 

C 

= 1 

.488 



= 1 

.25 


jr 

= 50 

.00 



= 58 

.59 

Kcal/M/h/°0 

T 

= 10 

.20 



Coefficient of linear expansion for concrete 
I = 11 X 10“^ per °C 

Mind Pressure; 

Basic V'/ind pressui'e 

100 kg/m^ from 0 m to 24 m 

115 kg/m^ from 24 m to 56 m 

p 

122 kg/m from 56 m to 80 m 
Shape factoi* = 0.7 

Earthquake 1 3 j 

Damping factor = P.02 
Factor =1.0 

Importance Factor I = 1 .5 
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(5) laterial Properties; 

Cube strength oP concrete(f ) 

Yi '.'Id stress of mild stecil (f ) 

sy' 

Unit weight of concrete 
Unit weight o.f fire-brick lining 
Young's modulus of mild steel 
Young's modulus of concrete 
Modular ratio 

(6) Material Costs; 

Cost of concrete shell ((\^) 

Katio of unit cost of steel and 
concrete (S) 

(?) Permissible Stresses; 

For H 250 concrete and mild steel as reinforcement 
the permissible stresses for warious load conditions 
have boon computed and are given in table 2 below. 

TABLE 2 

PERMISSIBLE STRESSES 


Load conditions 

Permissible 

stresses(kg/cm2 ) 

Concrete 

Steel 

JJ 1j , 4* * Tj , 

95.0 

1480 

D.L. + E.Q., 

100.0 

1560 

f . Ii . + Temp . 

82.5 

1430 

D.L. + W,L. + Temp. 

125.0 

' 1690 

D.L. + E Q. + Temp. 

125.0 

' 11690^^; V 


= 250 kg/cin^ 
=2600 fcg/cm^ 

=2400 kg/m^ 

=2100 kg/m^ 
=2.0x10^ kg/cm^ 
=181818.18 kg/cm^ 
= 11 

= Rs. 180 per m^ 

= 80 



• 5 io tij on Oompusor Proi^am: 


'hu caimney lias been divided inco 10 stages. In 
.gj:n. ral the stage boundaries are taken at corbel points, 
ii computer pro gcara in PGR IHkd is developed to solve the 
proseiiu problem, and the same is reproduced in Appendix A. 
Ihe computer program consists of one main program and ten 
subroutines. The main program initiates the program 
execution and prints the values of optimiun design variables 
at each decomposed stage of the chimney ^ The subroutines 
are briefly explained below: 


SUMi^ 


NEW CON : 


QUAD 


GRAD 


HESS I : 


This subroutine performs the constrained 
minimisation using extended interior penalty 
function method. 

This se.^gment is called by SUMT and performs 
the unconstrained minimization of the penalty 
function using Modified Nev^ton’s method. 

This modulo is called by Newton for performing 
one dimensional minimization. It uses 
Quadratic interpolation method. 

In this the gradient of the penalty function 
is evaluated. It is called by NE17T0N. ■ 

In this subroutine the search direction for 
linear minimization is generated after evaluating 



--"iiLYS 


AifECj 


Mil Vim 


PENAL n 


Ai'TEUr 


lessian matrix. It is called by NEW TON. 

Elis ae;5ment does the analysis of the 
chimney and calculates the values of tne 
constraints, 

Ihis subroutine eliminates the no^gacivc 
design variables before any analysis is 
performed, 

This modulus finds the inverse of the Hessian 
matrix and is called by HESS I. 

This subroutine evaluates the penalty term of 
the penalty function. 

It gives the value of p (Neutral axis of the , 
chimney) by solving the transcendental 
eq_ nation (2.21). 


4 » 4 Pertinent Values of Parameters used in the Computer 
Pro 'gramme ; 

fhe initial value of transition parameter (g^) is 
taken as 0.1 and the exponent (p) of the -penalty parameter (r) 
is t-'ikeii as 0.5. The initia.! value of r is chosen such 
that the penalty term equals the value of the objective 
function. After each unconstrained minimization r is 
reduced bv a factor of 10, 



0x13 dimons ion'll searches are conducted by using 
a p’!r'ibolic approximation to the penalty function along a 
scare :i direction. Ihis appro :iira':.t ion is based on the 
three bes c design poinrs previously determined when 
any of the three criteria are satisfied: (1) 20 steps 
have been taken in the search direction : (2) the 
expected improvement in the object function P(X,r) based 
on the parabolic approximation is less than 0,05 percent; 
or ( 3 ) the step size is reduced co less than 0.01 times 
the starting step size. The starting step size is taken as 
0,1 in each one dimensional search. 

C onvergenoe Criteria: 

In the presenc problem the iterative search 
process is Germinated according to the following convergence 
criteria. 

(i) j'or each unidirectional search 

P ’ ' ( X*, r)-P' (X* r) 

P’ (X*,r) 

ISspP j < 0.001 

(ii) for unconstrained minimization 

P ' ’ (X, r)- ^ ^ 

I P» (1, r ) 


< 0.0001 


< 0.005 



( i i i ) 


•or final completion of the search 

1 X : ' - 'x M < 0 , 01 

1 ^ 


F ' ’ (X) 


f 


(I) 


? > 


(X) 


< 0,0001 


In the above criteria the notations ' ' and ' 
refer to two consecutive points of comparison. 

4 . 5 Results ; 


Ihe results are obtained for each decomposed 
optimisation problem and they are compared with those 
in rtef, j2, 11 j and a,re tabulated in fable 5* Only 
mean diameter and wall thickness of th'S chimney are 
tabulated in fable 3 for each case, since the percentage 
of reinforcement at each cross section attains the 
'iiinimum peririissible value of 0.5 p-ercent a'c each optimum 
in all the thre'C cases. 



COMPARISON OP OPTIMUM DSSIG-I WITH THOSE REPORTED IN REP. l2,11 
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4 . 6 Bisou33ions : 

It can be seen that at the optimum design, the 
thickness of the shell at all sections reaches the lower 
bound value and so does the value of the percentage of 
reinforcement at each cross section, The mean diameter 
of the chimney from the top upto the flue opening remains 
constant, and it increases to a larger diameter at the flue 
opening and remains same upto the base. This has been 

i 

consistently obtained for this structure by using other 
optimization procedures in Ref , 

i*.t all the cross sections, the active constraints 
are the side constraints on thickness and mean diameter. 

The most critical section, for the chimney considered, is 
obviously the section at the flue opening. The behaviour 
constraints, viz., compressive stress in concrete due to 
dead and lateral loads, and due to dead load, lateral load 
and temperatures, were observed to be critical at the flue 
opening. 
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The computation time required to solve the 
present problem is 40 seconds in D30 1090 system. It 
took: 5 minutes of CPU time Ref. {11 j and 10 minutes 
of CPU time Ref. jOj for the same problem on IBM 7044. 

-But it is believed that the present computer (DSC I090) 
is faster by a factor of 4 to 5 as compared to IBAI 7044. 
Had this present problem been done on IBM 7044, it would 

. s , 1 

have taken 200 secs (40 x 5), that is approximately 3^ 
minutes. So it can be concluded that the present 
technique to find the optimum solution is definitely 
superior to those used in Ref. 12,11 1 , for solving 
this class of problems. 
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to^,al TninlTiiim cost obtained In the present work as 
Vifc-.l]. as in 'iof, are almost the samo. 

Las rly but not the least, the penalty— function 
teclinique seens to be gainiinp increasing popularit:/ in 
the field of structural optimiaation, and the Tolume 
of literature on the practical applications in this 
field is incre£ising rapidly. flie modified ITewton's 
method for unconstrained minimisation in the, penalty/ 
function forraulation, is robust in the sense it works 
v/oll though the iiossian matrix is not positi^re 

d 0 f in i t G . And ;al so 'wi t J!l ttlO inclusion of extended 
penalty .function technique, the raethod is able to accept 
infeasible starting points. Lhis is a quite important 

t 

feature when dealing with complex design problems for 
which feasible starting points are not readiljr aya.ilable. 
Although tlio method is already used v/ith success in 
practical design work, there are great potentials 
for furtac'i' improvements. 

4 • 7 Suvhostions for Further Work: 

(i) lost of foundo-tion could be included in the 

objective function along with the present 
objective for the minimum cost design of 
li.C. chimney. 



(ii) jierodynfimic effects on the chmney should be 
considered for more realistic analysis of the 
structure. 

(iii) lost of brick lining as well as cost of shuttering 
could also be included in the objectiro function. 
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r;)»-,io.|/,:i!pl/,c;TFP,pP31 ,FP,S2,MAXFIT,STI-H’’’ ■ 

rAG’^,''>TTM.:4Ro ’ ' 

<"uP^n*4 ffPF'^A ■ 
rfj*'';-ir'i.< phi^fw 

?.'Fr'/A»C ' , ■ ■ 

^?TTm #0 

W'fT^aMTIWft : 

I¥PK fl 5 ,KTlf 4 ^ ■ , ', 

■SX^'NP. or CT^Lr* % T2,5X/'» 

jr fm CYCliE TN.E^^CH UNCO^iSTRSTVED iWl^IUATlQU 

.,>^^mop'fp,r^hFh'i 

»a » w ‘ 

,. jK(i)+opTP*sfn 

’lM.ro, J) PHTOr.DaPAT.rA 

iriMTiw.FO n 00 10 

PHTNrwsP»i,rA : 

riiFfK r.'i)D f”0MVrP'^E0C'^ FSTN^ PHI VAOUE 


T}- f I rpr (Of.n.pmwEWI / phTUT.DI ,f.T. 0 . 005 'S ^9 TO 100 
Pl.l.n 0=PM1 ^'K'-.' 
t (• rii-r ; 'v.^CT.-n nu tOO 
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■ - r- .!HF frt UMli Dli^^iiSinigAb SFAi’CH 

■ '>‘.\rir Ibl’ePk’HLI'TTu.^ TfiCaNIOUF 

MMAr'<(r'E,S,.S,M,ALFA,FALFA) 

• F- ; fui) ptRiBj.S(J),nESl(3),DAf.Cl) 

l/aTh'P,FPSl ,Kro2,MAXFTr,STb»-?T 

■ > . ^ Hj'i A i ' , 

^ j ' (•■■), ZAP LHwi /Pat 

i';> ijf'h.VA ■ ■ 

F f K A s. (i . . 

fzrT^Ju ' ' ■ ■ ■ . ' ■ 

A TU , rt , ■■ ■ : 

ri t i> !•’ . ' 

b.,(p.i/-,<,ro'’'Trti?! PH' sf n , . 

i ^ 

l\i ‘y 

ni-. 1 f )n.f i) 

r j'*' i f i.ut.. 

?>5 1 (r>lj ^ ■ ■ ' 

DJ / J=l ' ' ■ . ' ■ ■ ", . 

<*» I, j ™ ^ vi ^ I 1 . ' . ■ - 

rji'i I'lwbk; 

■ ■ ■ ■ ■ - 

r.A*U«fcOr«tDKS,AI,S, A J ■ : 



■,W fo 20 



„ A*'*' 

MW4f:ho TO 120 

'■ 

RWAT 

F i sUNFOF« ( Oh A , r; , •?, , PT l 
rrf fi .hi.KAj v,p T,j .jA 
r,‘i rii to 

Fii-'- 1 


I , . ~ I H 

T\- ! ' ' .s .’f '-J v;**. Ti.i ^ 40 

1 I i ■' , 1 . t' . j ^ 'f'ij 


'h.l i ' ! P) 


I i’ll ou 
FH*H 
GJ i’a hU 

Ptm'MG ThF OUADMATXr 

C2«FB*(C«A1 

C3*FC»CA-lj) 




Aot'A=v.;>»lC4+C'^ + C!>)/(.Cl+C2 + C3) 
P ■V!.h h-0 mofn ( UF.S ft4fS fhhFh) 
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' j S' f |s*C 1 ♦ fC--\ ) 

"" :» *d<'i'-*c2#r Mtr 
' : !t:* i,tc^ </uFf 
^ ' u’'nri i/f>r/r 

i M- » ')l'>^u(<”A + H*AI.FA + *it ■ ■ 

'-■'i 4^, u^H»Af<KA 

■ f (HAi,KA j-f'Ai.PA)/FAl.rA) 

ii i r~) Fn’f I " ■ ■ 

n^H. /'t,A,FA B PR r,FC,ALFA,FALFA,HALFA,FRR 

M '4'/., I UF 17,4 j 

( K O'rr^nT.KftxFiT) an tq i4o 

!{•' » iFhK.C.p.tFSI ) .AN0,(DFAijFA.bT.FPS2) ) GO TO 140 
* Al.Fm) Co TH so 

i r ( r lFA) go xn too 

AsAt.Frt 
FA- i-'Aj>‘A 
nj f'G -s^i) 

1 1 ( FB.f,T.i''AuFn) 

( - A I, FA 
FC=FALFA 
CiJ XO ho 
.kMB 
fAssrB 


GO xn 90 






iikIfirA ' . ■ ' ■ 

GO xn ho > ■ . ‘ ■ 

Typt; \ 

FllPi.ATf >, (Oa , 'F,il>.CTin!4- unpS MOT DFGHFASh’ AGUMG THIS niRFfTiaW') 
AiJ'O.st . 

r'Al, I ... 

r. ) I ; ) I 0 . 

Sr.;, i,r"A __ . ■ ■ ' ■ 

•■'i' ' M ' I /, P^a', 1 lb at AbFASrAH=%F1 0,4) 

" ,1 ■'■!: WO , T KT'r , ,.PK’f..A ■ . . 

F /'• A i , IFa , 'SnABFi- OF OuAUR/vTir fits; * ,17, /lOX, 'TOTAL NO, OF 

1 ”''v I’jf'! .s ; ' , n/) 

tS; iO ^> 0 ’) . : 

roo ^V| f / 1 ox , 'cam tini' trap in ftftefn jumps'/) 

Pb'i UHS - ' ■■■ . . ■ ■ ■ 

FOO ■ ■ . 



”JS‘ih; i-OH niiVfcf.JPiiFi-lT 

’"iKh bRMr‘i!>hS,nij;,Gaj 

ra'.r{1u),!)Gan,1j,uGB(iO),ijGF(lOJ 
‘ B j ,G(U^ J ' 

' ( '-‘I ! I )f'* f'uNG f 1 0 J , UKS f J ) 

■ }'->• ii./cC'HG} /,irut!,fiVAK 

/ivSl AGF,MrtM,NSy . 

"..r-'- .'(../sFTHurVF'iyHf’,Hi'yj--L,NlYpF 

M /COiv 5 

^ . , ' ■ . ■ , 

,cunkA, am,Coc 

:.A'‘Mnri/cnr,5'//C1 ,C2,Ht ,S!GhA,THAF 

'Kj S J T < , •jPu*'* , ■ ' . 

j=.),U ■■ . ■ 

VoUt, [ ,!ru,0 ' , . 

'<J Mi jsi,r-:,'AK 
!Fn/F.‘.,(vU.F.y,t).U) r,U TO 6 
IFC J.FJ.NVAn) HjsO.lJUFSf J) 

[FCa,i;;0.?'»vAF) GO i>n ? 

^ ■ ■ 

50^ W"7 " 



^ fctf'ii * * If %.* 

)E«t J)KOKiJf 
:Mih hnAhY(DF.3) 
)D' 15 li = l,«rijr! 
)Gi'UIt l=Cut'pf 1 1 ) 

!''f ‘ ■ 

.1) MM 


>,,i ;n 




MV. 


( * , J ( I i f i )/('/, U^dJ): 


j 


I # 

, I 
t I 


’ ^ , 


% ri ! ' * ( 

)d 7’^ 1 = 

M I * J J 5 

rKf . Vvr-i- GO tu b7 

t I- (./I GO Tu 

Ir t .< '»’.7 hMO Cu in 

'>.Ji-Gu '^•VkAo, 0/( i7(rtP«=*2) 

Gu rn bvi . .. . 

I Ff J'M.'SiJj .GT.ThAF"* Go iO 
S U « S u I'H f 2 . U * C n N S C J J / i ' P A P - 


r ^ i i I -.G' v'J 1 

= I , " M' t<, 

, ' col, 

< r. 


55 


53 


AM))*UCAJ,J)/tnAP**2 
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' ‘ ^ UJ) .GT,TkAP) Gij i’O 55) . 

‘U.M/TPAP-I 

c .;,0«uT4-f2 + :^.U*UT-J .0)+uGt.J,Tj/rRAP + ’J ^2 

‘ 1'^ I f ‘it: ■ ■ ■ 

I ^ ■■ ■ 

' / i r 1 ■ ■ . 

i) ' ■ . •■ 

r >k HI pOM tL*f;.5Ci!!)'>i"UC*(1.0 + COSRA‘^yFS3)) 

^ /) = (pi 1 )«i'ruc*( t ,u+cnsRA*DES3 j ) ■ 

r’n-'H 1 )*CJS8A%-CUC 

! ) . I 0 W S /, " I » >v V A P 
G.,( l,'i=r A if 

C Uli t j. ■ . 

Rt'^'lJHi • 

i*:.ui • ■ ■ ■ , 



■’ll f I ! n -: 

~ 'T i ft 


n‘H uk 55 ; j . a -^ ^EVEhuPmnf 


V.> ; ui-i ljKf> (3,1), isAf! f J . 3 ) , f I 0 , 3 ) , G 
10),uPSf J) ,0(3) ,Ar;CX) , A 


J' 
* 1 , 1 " 


I / .. rtj *' ^ ^ AH 

* • 1 ‘ Hijll / mfr ^ f .,^ 1 1 f eh , cttrp K 

Hi , //, li |;| i ^/ pfl ||,^ 

n ■ 

ifi 


f 3) , 
AMGf i 


^.•,/C':hF.i;/(’1 ,ri!,Hl ,5;inAA,TK?vp 
f ».,/ r ^ 5 .|.<*,, v ./ p 1 jCOFRA , AM,CnC 
laj , K'VAP 
.''- J ’‘i ( ,uV ,\F 

S 4 ' ’ I t , U • * , u 


I*;-; 
-'' j 

t’ I 
0 j 
S l» 
0.1 
It' 


; f f 
r f 


r , 


.t)=(U<> 

Or, 

’-! , A ' VAk 
■^ - ! , 0 


1 = 0. ,, 

( M'l ypF . F 0 , 2 ) 


f;o 


, , TO 14 

fFf :<TyPi;;,rjO. n GO ru i5 
If <C0&.SC T) ,i.F,TR.»t») GU TO tb 

2 . 0 / r CONS t T J *4 3 ) ♦ { OG f I , J ) * 0 r, ( I 

"■'* ■?(>. ' ^ ' ■ 

.00 TO t 7 

I , a )* or,r r , K)/ (consct )** 3 ) 



K1 ) 


iST 

..’lo*’ 

ii?(TB 

|f)Wt2,0»nG( 1 ,K)4riG(,l , j3/tTOAP4‘#3) 


i » r ) GCi ,!!()»( 6 , o » truMS ( n / T*up 


‘I* 


n 

! I f' 


0 *' , 
. 1 . 
»S| 

1 1 , 

* i ' 

• I \ 

( 3, 
i ’* 
iO 

!( ! . 


^)•z 


f .J 


) 


"■Nf j 1 

; j =, r + r . lS ^<■^^‘!■• c ; s 3 *coc 

*.r .-FHA + r rjSC > J *COC 
^ f f'uP.- i j^'jOC 

1 i=.s^ r 1 , <;) 

ij ~. sA / n , i ) 

■>)=.■, A 
1 = 1 . ' • .O’ K 
. I = 1 , 1 

i , V ) ) + ,'■: nM(T , J ) 


. 1 , 1 , ! jtJMP ' O.f i , J ) 


TYM.: !v’. , ( I Of.'? ( T ,.,T) , jsi ,;jVAH) ,1 = 1 ,i<iVAR) 
FilH.iAT ( / H>X , ".if 531/. i -lATRlX; ‘ ,/3(iOX, 3 Cei 4 
HT’SF-lA-i MATRIX TT bDGT FORM 

no HO \^\.osrhH ■ . 

A,MG(Tj.sy,i) ' . 

AbtDaO.O ■ • ■ 

rH . f .)« o.o ■ 

CDlfni'rof ■ 

, t), 

AMGCI )«liKS(l,2) ■■ ■ 

ALf n ». AMC ; C 1 J/|>f ■ ■■ '■ ■■■ 

fX ?)' e 5 ' fS ( 2 , 2 )- Aifl 3 * AMGtl ) , . ^ / 

‘'iX(‘U7) = ‘lF5f 1 , .3) ^ ■ ; . 


A( 3 ) 

» , 'n 3 i , 5 ( i 1 


-1,0) 


.a,5X)/)/) 
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' N 1 1 ~ALf n»AMG(2j 

j/i)( n 
i '>=^w-G(^)/i.f^) 

' 'J l=t ,i'.V/\P 

' ' ■ , 

>'> i , 

’* j'.' n >nr ■ . ■ ■ 

■’ ) 4 ^ 4 . us I , AVAP 


n-’Ci.r, , ?),!/!', i-.O) GO TO 2U5 

!'-j ; ? ■ ■ 

rVi^l, r ) , j =1 , jTl 

Frisj ?■. f/,sy, '■'.M T,>iVMf‘>X,P14.4l/) 

(■t/IMP.. 

'■, . 1 ,; T’-n M'-rr/ 

fv.) 2 ('7 i = i,u'’i'.W 

Jl- ( Jt i ),hi .u,Ol A(T jsl.tj 

TCJ)=/\f 3 ) ■ ■ . 

T^2l«AC2?-Ai.f J)^TO) 

TCI }»A(l)-ALn)»T(?)-At4(2)»TC3'> 


TH 240 


"rtoK- ' " 

TlfPK 500, (F) r j , 1=1 , i) 

FiJW'T.'rf /sx, *>'"n DIRK: ' ,3t5X,E1 4.13)/) 


^ nil I’H 

7 5 0 

1 . 1 !'|1 

1 = 1,. VAP 

n r jr:. 

, i 1 I ' ; 

Mil 

•■'!• 

'1 } ^' > ^ 


r 
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KUNCTxnh, Of^KnF^tDER.N^S, ft) 

,0ES1 C3) .S(3) 

rov^PN/rowSa^CI ,C2,K1 ,SI^MftpRftP 

rowMnr./AREft9/cnws 

ro«^Mnrj/AREft/nESl 

roMMHM fyPEVA 

RO in Tsi.w 

nESlCDsnfesCD+A^SCl) 

type^^oesi 

TYPE*, ft 

00 !*> TsI.M 

TFCDFSt (T) ,GT,0,0) GO TO 15 

GO TO 17 

roWTTNPE 

00 TO 18 

GAOL ANEOCOeSl,S,EK) 
ftsEK 

no IP J=1,N 

OESl CJ)=nESCJ)+EK»SCJ) 


nONTTNOE 

CAl-G Aft'Ar.YCDESl ) 
CALL PFNALTpESl. 
OMEDENbORJFN(DES1 
TYPE^.ONEDFN 
PETURW 


FCONS) 

^ND+FCOivJS 


ENT) 


function DFUN(DES,N..5, A) 
niMENSTON 0E5-C3) .sr I) ,G( 1) jDESl C -' 
ra»<MnN/CPN*52/Ct ,C2,Rt ,5IGMA,TRftP 
roMMON NFEVA 
no 3 1 = 1 , N 

OESl fI)BnE 5 (T)+SfI)*A 


rONTTNUE 

gat.l GPftncrtKsi , do,g) 


F-o.n 

no 5 1=1, N 

' 5 '=F+SCT)*G(:i) 


HFUNsF 

PETUPN 


function 0RJFN(DES,N) 
niMENSTUN HERO) 

COWMON NFEVA ... 

rOFMON/P ARAM/PI ,COSRA , AM ,COC 
COmMOM/ET/AH 
n p SE n F ^ f ) 

OBJFN = (PT*nERC 1 )*DFSC2)’»=(1,0+C0SRA*DES3))*C0C 

NfFVft=NFFVA+l 

RETURN? END 
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v.';ic- 0^ CHi«>'fK'y- 

■’.riijf (i>fs) 

■: -‘jJ' ».'KK.4'1 Of CC‘-Ci':fc' re 
■■.<- !/.!/»', tTFK tjF C!.l,'.wxex at TDi-* . 


Ti.H uiA-’ETKt? Of Zhl AT fiV 

'!■',- s ru-.KA .:a/,.iFTF:K OF Zi^ ■'■i'-iW aT UF 
'-■;‘0'!iKF i-i A.''-FTt:F UF Li;Mi4G A'f J JF 
ii< i.lt i >,-'.k.H OIAp>Fli.'tF OF LlivlfiG AT SFCtiO;'? 


J lUf-.H (!)) , A ii'i'J' i'f ' u 


SbCi it'A CO ! ! A''> M * ;*i 4 V H /’I J t'.' 




( -U'. *'. r. * I j' 


■Ji oJa of LT.’a-VG .AC Ti)F 

wKi'iF'rHr ink nr i.jwi?-'c: nio 

.CO'is ihi.OM'l (if- CDiK^Ki F iGiFtO 
sSf,.K!G}ir (iF LVuluC, 

;jFs. FACT lU' 

Mf. -iV'FAA jOAM-.TFR of aT HDTi\J.‘i t SFC F H 4. ij 

FiA COlMA'Fi ShEiuF AT ThF AHOVsl StCliF 
l'{A OF iHtOU Al I.iri AiiOjK GOC 

OJ;«-FAA UlliWFTKP OF aIH SPACF Xu aFI 

AU'tiN$|00 £?ESC3},rti?^vS(U') 

CO.-iiV*0,*/AHeA9/CO<^S 

CaOA0F/At?f:A4/w»A^i0M,ALP,rtU,ftC, aCD 
^i:i;'ir,0A/AHF.A7/0{)BG(tl), JTH(1 i >,SAb<.%( F) i 
Caf«Ow/A*^|!A3FTOTl.,tOI^,AlHGftF 


1 3 



:AS/p.,PMa4KTA r'e.tna 


I.) /' r H 


/tH 

lAN/PI ,COSHR., A-’-uCOC 

»,THIT ,i)(Ui>F,HSl,:;'*iVU,22S,a, l «,.•< .^78,':3u,U/ 
(» A P f 4 > ) t \ ^ C 0 / F w 0 ^ T * i (f 0 » 1 Ji S ^ 2 4 s) u » 0 if 2 i U u w V / 

SP , OifH , r AX , FOOT , Ho/0 , 1 0 , y . U 7 . i o . !j , C , 5 / 
CZfZ ,M- 1 , AK/, 'i-'-S/l ,^»HM,J ,2A,5j,a,5H.S4, Iv,;,/ 

■ r , ■ '■» r F" t , 1 / 1 5 4 1 o , 1 d * 2 F t) , i 1 F " fo / 

■ - , f M, ., F/v ,i’/i , i .0/ 

'" , i-’ ' A / '/‘lO , V* f ?ns,»0/ 


/ 3 f, . 


rtI».0,U.U/ 


fA.M G 
i ) 


J l'( I I 0 
roTOsHETG'iT 


; (/'i .,s, * .OF PiiFSKfVl SKCIIOII XS Af .A Dr^FTAiJCF OF %F4,1 
I X , * '-e.riOvS F, '!!'•> rpHj BO (’TLli''i V) 

f I' ' c J A r! L’ I ■ ■ ■ ■ ^ ' ' ' 


i i - , . 

AlACsrOUATC .SlAGFJ 
GiAHflGvT Ur.TXGtt VAPJABuFS ..A 

y|,SCI)aiil‘fAXNC*0.4 ^ 

ifd-iSTAiJfi.le.S) l)FSC2JsO,it> ■ , , . A- 

IFC (W8fAGe.a|,6)*Aiin,CNSTAa€.GF, ; ) ) OFS(24=U.2 p 

If (.4§t^x;E,G'E*85 275 , „; , 

'w;siix»off05'A 


oqDir'Tgl) ArtAtysiS^ PROGRA*^. 


Jeiil 


# 


■ es 

/ 


' l '»H'l "Tu’f ij ■ • 

1 ) GO ro 

t ' . 1 SOu 

■ivl.4STAGP>U ' ' ■ ■ ' 

1 4 - !'»» ) ■ - ■ ' ' 

I : ^ . ; ' ; . . ; ■ ■ , . 

■: ' ■ : . ■ 

(j*'ALlfv . ■ ■ ■ ■ 

II' ■ ■ ■ ' ' ■ ■ ' ■ ' ^ ^ ^ " 

DBOsad+Tii r . ■ ' ■ " ^ . 

lii' 

t)l' li. , . ■• , ■ 

USsysi.u + AiRGAp ■■ . 

i,5!’l/f»i^(4l*a«'2,0»TjuJ'M • ■ V. ■ . ■ ' . : ■ ■ ■' 

A-!AiiCi'»t}6(iii}Bt*+DTo*DTa4Dan*nTO-jBi'* jri t-Dti »im- ).M4.)'n 
y^C,(i"4l*isAMASC'i'*F'CU»AH*PI/12.u , , 

A'''iASijT»0aLO*DBLO4DLO*OuG + UMin*ULO-OBL i i 4i/t,J ..u, ( 

B * •'CQ'itf^WiawT 

+ «C(iNT ■ ■: ' ■ 

AaiS.aT.fc) GO TO 42 , . ■ 

||MASL‘l+l»BLO*GGO»OHGI»Oul4 2.0 + lOi(U*BOO»i.'Li H)ht J 
itA*'it)f''-l*AH*0.25/A'dAGb’V 

MASC 1 + 2, 0+ ( uTj + ii f O + (l)B;j*irf!>-i)h i 4i) ri 5 

sa-m-./a AS{:’r+a,2r-.^Mh ' ■ 

Cfi K:h i 'Iff I ’ sr.c’i'T.v.., . . 

t'-= ti S( I j ' ■ , ' 

f } ' ■ . 


t ' , . 1 1 , 


T 


: u. ■ 

( 

A.-i' 

A'JU 


K 


■ -,J' i: » 

,,, A i i, 

) ' j i 4 Z ,t 0 - : 

\ 1 ‘ ’ *■ i ' I i 4 :• 

li' * 4 '! > Y 

^;>L''4 , , 4 V i 

■ ^ 1 - I :/ . •. ♦ ( ; 

i . I ! i,.4 M . iJ. > 1 +'.f 
= '<"ljrJ-4A. I' 

I 4 OiCtiAHTf 


( hP . ■ ■ 

I i.T .'> ■ ^ - • 

i 442 ~ PA:U 4 + 2 -i) 0 I 4 + 2 + DBBn 4 'DHl)-D 8 BJ. + l}Hl 


♦ 42 + 40 HbI-DBiy, i 

♦>' 2 b / -I 2 , 0 \ 

'Y« 4 i 4 42 ) + 0B04iMrttJ-i)BI*UBHi 

* 42 l+DHBL.j 40 rtb(’-UHBnl 40 BL!: : 

?H/,\ ., uSL’B 
2 B 2 A Y A hj 0 « 

0 T b ) + A G I i'i T* ( fi i- rt A rl T 4 TOT ) 


/i;..iJ + ACU!^B4HCBAHH4Wia>»3*‘iL!:UHB + 4* J4D ,e>4 cot'll i;HT+T., I fl. ) 
l)iYft«WCaN'i’ 4 hl.I(’i'J 44 »bJ‘'ief'WC 0 ' 4 'rt + w.J . ■ ■ ; ; 

rieAf<*AbU/Dli,w, _ . . ■ ■ ■ ■ 
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GO i'"i VL- 
,Q '('J l(') 


' ■' ■' . i*) t,(f TO ‘>uiu 

■ ;• 'O V'O- vO l r )-.SAByGl 1 H. 3 ) 

' ' ‘ I t> 1, 4 y 

i i f -3 w 5 tj ■ ^ ■ 

• ! : 

> f s'irt'i -“fC^Ti. 

'iO’s'ff / I , ' ■ ■ 

3 t v' , 6* i AL-i/TjO + ♦ J , H fy , 4-*^ ! A). 

O'./ ‘t’fj 2UO ■ . ' 

<M j ,5o,0),3’-.0. (TOtu.GS.24,yj ) 

It 1 ( i 0 1 1 » ^ Ci I # 2 4 ^ 3 ^ A 0 * C f * f G * 0 c, « # 0 } ) . 

'I' 1 K 4 '5 

c "EKf-Of Id CAi.CdLATJ 4^; 

EH i # t .■■.■ ■/ 

iiu TU i ^1 u 

PHEls|2t^.i) , , ■■ ' ■ . ■■, 2', . „ .' -■ ■■ .or,., ■■■ 

PRKaaJl^.O '■■ ^ ■■ ■ . ' . ■ 

n'dmfiis*h4,Q . ' . ■ , O' ■ 0 , ■ ^ 

y!?Hf:i*0.5»PREI*(&.O+QDBO( i,l+DrH(3) J^SP + Hi . 
'/»PHt2»a.5*PP,£2»(iiul^OC3)+Ui‘H(3)+UBO)^=Sr’*d;. : 

? W*A‘PK6i ^ { HI / 3 , 0 ♦ ( 1 1 0 . ')4 iiOfiO ( 3 ) + JTti ( 3 ) / / ( 5 , ^ + 1 ly toi ( J ) ; | ■■' ( y i j j 

^Mp4|«W,PRe2*i‘{H2/3.0»{'a,a*(UDBa{J3 + OrH( 3)3 »-OBO)/(eou + i.iL'OOH U3 

♦ AMUM2 ■ ' ' 

OfO- 2Cy ■ 


Ort/lOA, IG mT iih CHl.'^.vKY ABOVE THE PREBfv.T BE 


r. 


I. Ci i ' I t‘ J *'. f ? I 


f 4 , < 

H '‘>Z > / . . 


."sr. nj*'i .J Kir n .1 ■■! -ti K. i hc^kil inn^ r irv Cl, ' ^ I vj c, s, i • -- - -- 

'EC,', / lox , * 1('1‘ AT That. tSEGT I J'l I S ' , E I / , 2 , 1 a , KG« -iH ( ' / ) 


• • ~ ( y .‘J . . . ■ 

' 0 = ’ I 0 . . '0 ' ■ ■ ' 
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